Baryon and Antibaryon production at RHIC energies in the Dual Parton Model 
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We compute the mid-rapidity densities of pions, kaons, baryons and antibaryons in Au-Au col- 
lisions at yfs — 130 GeV in the Dual Parton Model supplemented with final state interactions, 
and we present a comparison with available data. 



1 Description of the model 

The rapidity density of a given type of hadron h produced in AA collisions at fixed impact 
parameter, is given bytH'H 



dN AA ^ h . - „.,r_.,i i 



dy 



(y, b) = n A (b) [JV« ( "* (y) + N*$* (y) + (2k - 2) JV*"Jj (y) 

+ (n(b)-n A (b))2kN^(y) . (1) 



Here n(b) is the average number of binary collisions and n A (b) is the average number of par- 
ticipant pairs at fixed impact parameter b. P and T denote projectile and target nuclei, k is 
the average number of inelastic collisions in pp and p(b) = ku(b) with v(b) = n(b)/n A (b) the 
average total number of collisions suffered by each nucleon. At y/s = 130 GeV we have k = 2!. 
The Nh^g,) (y) in eq. (|lj) are the rapidity distributions of hadron h in each individual string. In 
DPM they are given by convolutions of momentum distribution and fragmentation functions @. 

It was shown inll that eq. (|l]), supplemented with shadowing corrections, leads to values of 
charged multiplicities at mid-rapidities as a function of centrality in agreement with data, both 
at SPS and RHIC. Here we use the same shadowing corrections as in ref. i 
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Net Baryon Production (AB = B — B): In order to reproduce the stopping observed in Pb 
Pb collisions at SPS, we introduce a new mechanism (diquark breaking) based on the transfer 
in rapidity of the baryon junction ou. In an AA collision, this component eives the following 
rapidity distribution of the two net baryons in a single NN collision of AA u 



dN DB 



dy 



(y) = c m [4 /2 (i - ^+r (6) - 3/2 + ^ /2 (i - z^-y*} (2) 
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where Z± = ex.p(±y — y max ) and u(b) = n(b)/riA,(b). C v (p) is determined from the normalization 
to two at each b. 

In order to get the relative densities of each baryon and antibaryon species we use simple 
quark counting rules a. We denote the strangeness suppression factor by S/L (with 2L + S = 1). 
If the baryon is made out of three sea quarks (which is the case of pair production) the relative 
weights are I 3 = 4L 3 : 4L 3 : 12L 2 S : 3LS 2 : 3LS 2 : S 3 for p, n, A+S, 5°, S~ and Q, respectively. 
The various coefficients of ^3 are obtained from the power expansion of (2L+S 1 ) 3 . For net baryon 
production there are two possibilities: one is to use 13 which corresponds to the transfer of SJ 
without quarks. The other possibility corresponds to the trasfer of the baryon junction plus one 
valence quark. In this case the relevant weights are given by I2, i.e. from the various terms in 
the expansion of (2L + S) 2 . This second possibility is favored by data and it is the considered in 
this work. In order to take into account the decay of S*(1385) into An, we redefine the relative 
rate of A's and S's using the empirical rule A = 0.6(S + + S~) - keeping, of course, the total 
yield of A's plus S's unchanged. In this way the normalization constants of all baryon species 
in pair production are determined from one of them. This constant, together with the relative 
normalization of K and n, are determined from the data for very peripheral collisions. In the 
calculations we use S = 0.1 (S/L = 0.22). 

Final State Interactions: The hadronic densities obtained above will be used as initial con- 
ditions in the gain and loss differential equations which govern final state interactions. 

dT u k 

The first term in the r.h.s. of (|3|) describes the production (gain) of particles of type i resulting 
from the interaction of particles k and I. The second term describes the loss of particles of type 
i due to its interaction with particles of type k. In eq. (||) p\ = dNi/dyd 2 s(y, b) are the particles 
yields per unit rapidity and per unit of transverse area, at fixed impact parameter, a^i are the 
corresponding cross-sections averaged over the momentum distribution of the colliding particles. 
The channels that have been taken into account in our calculations are 



vrN ^ if A(S) , vrA(S) KE , ttE KQ (4) 

Of all possible charge combinations in these reactions, we have only kept those involving the 
annihilation of a light q-q pair and production of an s-s in the s-channel. We have also taken 
into account the strangeness exchange reactions 

ttA(S) i± KN , vrS ^ KA(H) , nVL ^ KE . (5) 
as well as the channels corresponding to (|j) and @ for antiparticles. 



2 Results 



The calculations have been performed in the interval —0.35 < y* < 0.35. In Fig. la-Id we show 
the rapidity densities of B, B and B — B versus h~ = dN~ jdr\ = (l/1.17)dN/dy and compare 




Figure 1: (a) Calculated values of rapidity densities of p (solid line), p (dashed line), and p — p (dotted line) at 
mid rapidities, \y*\ < 0.35, are plotted as a function oLdN~ /dr), and compared with PHENIX datau ; (b) same 
for A and A compared to preliminary STAR datau ; (c) same for H~ and 3 + ; (d) same for 57 and . 

them with available dataBi'i. We see that, in first approximation, p, p, A and A scale with 
h~ . Quantitatively, there is a slight decrease with centrality of p/h~ and p/hr ratios, a slight 
increase of A/h~ and A/h~ and a much larger increase for H (H)//i _ and f2 In the DPM 

(before final state interaction) the rapidity density of charged particle per participant increases 
with centrality. This increase is larger for low centralities a. This has an important effect on 
both the size and the pattern of strangeness enhancement in our results. It explains why the 
departure from a linear increase of H's and fi's (concave shape) seen in Figs, lc and Id is also 
more pronounced for lower centralities. 

The ratios B/B have a mild decrease with centrality of about 15 % for all baryon species - 
which is also seen in the dataK Our values for N ch /N^ ax = 1/2 are : p/p = 0.69, A/A = 0.72, 
S/H = 0.79, Q/fi, = 0.83 to be compared with the measured values! : 

p/p = 0.63 ± 0.02 ± 0.06 , A/A = 0.73 ± 0.03 , S/H = 0.83 ± 0.03 ± 0.05 . 

In Fig. 2a and 2b we plot the yields of B and B per participant normalized to the same 
ratio for peripheral collisions versus n par t- The enhancement of B and B increases with the 
number of strange quarks in the baryon. This increase is comparable to the one found at SPS 
between pA and central PbPb collisions - somewhat larger for antibaryons. Before final state 
interactions, all ratios K/h~, B/h~ and B/h~ decrease slightly with increasing centrality. This 
effect is rather marginal at RHIC energies and mid-rapidities. The final state interactions (Q), 
H lead to a gain of strange particle yields. The reason for this is the following. In the first 
direct reaction (||) we have P-k > Pk, Pn > PA, P-kPn 3> PkPA- The same is true for all direct 
reaction (||). In view of that, the effect of the inverse reactions (Q) is small. On the contrary, 
in all reactions (||), the product of densities in the initial and final state are comparable and 



100 ISO 200 250 300 350 400 



50 100 150 200 250 300 350 400 



Figure 2: Calculated values of the ratios B/n par t (a) and B/n par t (b), normalized to the same ratio for peripheral 

collisions (n part = 18), are plotted as a function of n par t- 



the direct and inverse reactions tend to compensate with each other. Baryons with the largest 
strange quark content, which find themselves at the end of the chain of direct reactions (||) and 
have the smallest yield before final state interaction, have the largest enhancement. Moreover, 
the gain in the yield of strange baryons is larger than the one of antibaryons since ps > Pa- 
Furthermore, the enhancement of all baryon species increases with centrality, since the gain, 
resulting from the first term in eq. @, contains a product of densities and thus, increases 
quadratically with increasing centrality. 
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